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CHAPTER  I 
INTRODUCTION 

In  this  dissertation  we  consider  some  problems  of  inference  for 
accelerated  life  tests.  Accelerated  life  tests  are  used  to  obtain  fail- 
ure data  on  items  of  high  reliability  in  short  periods  of  time.  An  in- 
troduction to  the  topic  of  accelerated  life  testing  is  given  in  Chapter 
9 of  Mann,  Schafer  and  Singpurwalla  [11]. 

For  the  purposes  of  this  dissertation  we  consider  parametric  ac- 
celerated life  testing  procedures  which  assume  that  items  under  test 
have  a particular  failure  distribution.  We  also  assume  that  the  param- 
eters of  this  failure  distribution  change  with  the  stress  according  to 
a prescribed  functional  relationship.  That  is,  the  parameters  of  the 
failure  distribution  may  be  re-parameterized  in  terms  of  the  stress. 

The  re-parameterization  appropriate  for  a particular  life  test  is  de- 
rived through  consideration  of  the  failure  mechanism  of  the  test  items. 
It  is  understood  that  a re-parameterization  model  is  valid  over  a range 
of  environmental  stress  that  includes  use  conditions,  or  normal  stress, 
and  accelerated  stresses.  The  accelerated  stresses  are  large  enough  to 
induce  failure  in  a short  period  of  time.  A similar  situation  prevails 
in  certain  kinds  of  biomedical  experiments  wherein  the  parameters  of 
survival  distributions  are  re-parameterized  in  terms  of  concomitant 
variables. 

This  dissertation  is  divided  into  two  main  parts.  In  the  first 
part  we  consider  a non-Bayesian  approach  for  the  analysis  of  three  mod- 
els which  re-parameterize  the  scale  parameter  of  an  exponential  distri- 
bution. These  are  the  power  rule,  the  Arrhenius,  and  the  Eyring  models, 
respectively.  In  the  second  part  we  consider  a Bayesian  approach  for 
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the  analysis  of  accelerated  life  test  data  in  which  the  power  rule  is 
used  to  re-parameter ize  the  scale  parameter  of  an  exponential  distri- 
bution and  a Weibull  distribution. 
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CHAPTER  II 


NON-BAYESIAN  ANALYSIS:  LEAST 


SQUARES  ESTIMATION 


1 . Introduction 

In  this  chapter,  we  consider  least  squares  estimators  for  some 
accelerated  life  test  models  involving  an  exponential  distribution. 

The  models  considered  are  the  power  rule,  the  Arrhenius  model,  and  the 
Eyring  model,  and  may  be  written  in  the  general  form 

0^  = expCX^B) 

where  0^  is  the  scale  parameter  of  the  exponential  distribution  at 
stress  level  i , X.  is  a vector  of  stress  values  that  represent  the 
ith  stress  level  and  6 is  a vector  of  unknown  parameters.  The  gen- 
eral form  of  this  function  leads  to  a linear  regression  model  which 
expresses  the  relationship  between  log  CL+  and  B . Data  from  life 

tests  conducted  at  several  accelerated  stress  levels  may  then  be  used 
to  obtain  least  squares  estimators  of  B . 

Various  approaches  for  estimating  the  parameters  of  the  models 
considered  here  have  been  studied  for  different  life  distributions.  Maxi- 
mum likelihood  estimators  for  the  power  rule  are  presented  by  Nelson  [12], 
Singpurwalla  [19]  and  by  Singpurwalla  and  Al-Khayyal  [21],  The  maxi- 
mum likelihood  approach  was  also  used  by  Singpurwalla  [20]  for  the 
Arrhenius  model  with  an  exponential  life  distribution.  Graphical 


Note:  All  logarithms  are  to  the  base  e 
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methods  for  the  Arrhenius  model  with  a log  normal  life  distribution 
were  given  by  Nelson  [13].  In  another  paper.  Nelson  [14]  presented 
graphical  methods  for  the  power  rule  with  a two  parameter  Weibull  dis- 
tribution. Linear  regression  methods  for  the  Arrhenius  model  with  a 
log  normal  life  distribution  were  given  by  Nelson  and  Hahn  [16]. 

In  the  approach  we  consider  here,  we  define  linear  models  for 
the  re-parameterization  models  using  a log  gamma  error  distribution 
and  allow  the  use  of  censored  samples.  This  approach  follows  that  of 
Singpurwalla,  Castellino  and  Goldschen  [22].  They  obtained  least 
squares  estimates  of  the  Eyring  parameters  and  assumed  them  to  be 
approximately  normal  on  the  basis  of  a simulation  study.  In  this  chap- 
ter we  obtain  orthogonal  least  squares  estimators  for  the  parameters 
of  all  the  models  considered  and  establish  approximate  normality  for 
the  estimators  in  terms  of  their  exact  measures  of  skewness  and  kurtosis. 


The  methods  presented  in  this  chapter  depend  on  the  results  of 
life  tests  conducted  in  the  following  manner.  The  ith  life  test  is  con- 
ducted by  placing  n^  items. on  test  with  the  stress  level  fixed  at  the 
ith  stress.  The  sequence  of  stress  levels  at  which  tests  are  conducted 
and  the  assignment  of  test  items  to  stress  levels  should  be  randomized. 
This  is  done  to  help  Insure  independence  between  test  results  at  the 
various  stresses  and  to  avoid  possible  systematic  bias  entering  into 
the  test  results.  If  the  test  is  terminated  after  r.(<  n.)  failures 

l—i 

occurs,  a minimum  variance,  unbiased,  efficient  and  sufficient  estimator 
of  0^  is  given  by  Epstein  and  Sobel  [6]  as 


0. 

l 


r . 

i 

£ 

1 = 1 


t . . 


+ 


<Vri> 


r , 
i 


where  t^.  = jth  failure  time  (j=l,...,r^)  observed  at  the  ith  stress 
level.  The  distribution  of  9^  is  a gamma  with  scale  parameter 
(r^/CO  and  shape  parameter  r^  . 

Linear  models  of  observations  for  the  re-parameterization  models 
considered  here  may  then  be  formulated  in  terms  of  the  random  variabLe 
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log  0^  . Since  0^  has  a gamma  distribution,  log  0^  has  a log 

gamma  distribution  which  for  large  values  of  the  shape  parameter  is 
approximated  by  a normal  distribution  (see  Johnson  and  Kotz  18]  p.  196). 
Thus,  a linear  model  based  on  log  with  r^  sufficiently  large 

may  be  considered  to  have  an  approximately  normal  distribution.  The 
mean  and  variance  of  log  0 are 

E(log  0 J = log  0^  - log  ri  + iKr^)  (2.1.1) 

Var (log  0.)  = ij;(  )(r1)  (2.1.2) 

( ' ) 

where  !]»(•)  and  ij>  (•)  are  the  digamma  and  trigamma  functions, 
respectively.  Digamma  and  trigamma  functions  may  be  evaluated  using 
formulas  given  by  Abromowitz  and  Stegum  [1],  p.  258).  Thus,  when 
the  number  of  failures  r^  is  known,  the  variance  of  log  0^  is 

also  known.  Vie  may  Chen  define  a linear  model  based  on  log  0^  with 

the  observations  appropriately  weighted  by  the  known  variances.  The 
variances  and  covariances  of  least  squares  estimates  obtained  from 
such  a model  will  then  be  known. 

Since  the  variances  of  the  observations  are  known,  we  may 
construct  linear  models  that  result  in  orthogonal,  or  uncorrelated, 
least  squares  estimates.  A number  of  investigators,  using  different 
approaches,  have  determined  that  orthogonal  designs  for  linear  regres- 
sion models  have  desirable  and  optimal  properties  (see,  for  example, 
Ehrenfeld  [5],  Elfving  [6],  Himmelblau  [ 7 J , pp . 234-237  , Tocher  [23]. 
In  addition,  orthogonal  designs  greatly  simplify  the  process  of  com- 
puting least  squares  estimates  and  circumvent  problems  that  arise  if 
the  covariance  matrix  is  ill-conditioned.  Orthogonal  estimates  for 
the  accelerated  life  test  model  parameters  may  be  obtained  by  slightly 
modifying  the  form  of  the  model.  These  modifications  make  use  of  the 
known  variance  structure  of  the  corresponding  regression  model  and  do 
not  alter  the  basic  character  of  the  model.  The  least  squares  esti- 
mates for  the  parameters  of  the  models  considered  here  have  some  other 
convenient  properties:  it  is  possible  to  obtain  least  squares 
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estimates  that  are  best  linear  unbiased  and  unrorrelated  for  any 
sample  size,  i.e.,  any  number  of  stress  levels;  the  exact  variance 
and  higher  moments  of  the  estimates  may  be  obtained;  approximate 
normality  may  be  established  for  the  estimates  for  a small  number  of 
stress  levels,  provided  a moderate  number  of  failures  is  observed  at 
each  stress  level. 

In  the  following  sections  of  this  chapter  we  consider  the  power 
rule,  Arrhenius  and  Eyring  re-parameterization  models  in  detail.  In 
each  case  we  define  regression  models  that  result  in  orthogonal  least 
squares  estimates  for  the  model  parameters  and  develop  procedures  for 
inference  at  use  stress.  The  procedures  are  illustrated  using  some 
data . 


2 . Power  Rule  Model 

For  an  exponential  time  to  failure  the  power  rule  assumes  that 
the  scale  parameter  0 is  proportional  to  the  Pth  power  of  the  stress 
V . The  power  rule  for  0 at  stress  i is  written  as 

n.  = c/v.p  (2.2.1) 

1 1 


where  C and  P are  unknowns.  Life  tests  conducted  at  k distinct 
accelerated  stress  levels  V.  , i = l,...,k  , as  described  in  the  first 

l 

section  of  this  chapter,  provide  data  that  may  be  used  to  compute  the 
estimates  0^  , i = ],..., k . This  data  may  also  be  used  to  estimate 

C , P and  predict  0 , the  mean  life  at  use  stress. 

We  may  formulate  a linear  model  for  the  power  rule  that  results 
in  orthogonal  least  squares  estimates  by  slightly  modi fying  the  model 
of  Equation  (2.2.1).  The  modified  form  of  the  power  rule  is 

0 = C*/(V./V)P  (2.2.2) 


where 


V 


(2.2.3) 


■ 


- 7 
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with 

w^  = ^fr^)  * * = • (2.2.4) 

In  effect,  the  modification  is  a change  in  the  proportionality  constant 
C in  Equation  (2.2.1),  indicated  by  the  notation  C*  in  Equation 
(2.2.2).  The  value  of  P is  unaffected  by  the  modification. 

If  we  assume  a multiplicative  error  term  for  the  model  of 
Equation  (2.2.2),  a model  of  observations  for  the  power  rule  may  be 
written  as 

where 


X 

X.  = [1  , - log  v./V]  , i = 1,  . . . ,k 
6’  = [B0  , 81]  * [log  C*  , P] 
e'  = [log  ex  , . . . , log  ek] 

Z — [ , . . • , 2^] 

z±  = log  0i  + log  T±  - i[/(r1)  , i = 1 , . . . ,k 

and  E'  denotes  the  transpose  of  6 . From  Equations  (2.1.1),  (2.1.2) 
and  the  randomization  of  the  life  test 

E(e)  = 0 

E(e'e)  = W = diagfwj^, . . . ,wk] 

The  w^  are  given  by  Equation  (2.2.4)  and  W is  a k x k diagonal 
matrix  which  is  the  variance-covariance  matrix  of  e . 


■kw. 


3.  Arrhenius  Model 

The  Arrhenius  model  usually  has  applications  in  tests  where 
temperature  is  the  accelerating  stress  and  it  assumes  that  the  exponen- 
tial scale  parameter  may  be  re-parameterized  as 

A+B/T 

0.  = e 1 

i 

where  L refers  to  a particular  stress  level  and  A,  B are  unknown 
constants. 
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The  natural  logarithm  of  the  estimates  0^  may  be  used  to 


construct  a linear  model  for  the  Arrhenius  parameters  that  is  similar 
to  the  linear  model  defined  for  the  power  rule.  Orthogonal  estimates 
for  the  Arrhenius  parameters  may  be  obtained  by  slightly  modifying 
the  basic  form  of  the  Arrhenius  model.  The  modified  form  is 


0i  = e 


A*+B(T11-T) 


(2.3.1) 


where 


= v\1t±/ i,  wi1  * 


k 

I 

i=l 


(2.3.2) 


with  the  w^  defined  as  in  Equation  (2.2.3).  The  modification  is. 


in  effect,  a change  in  the  proportionality  constant  indicated  by  the 
notation  A*  . The  value  of  the  parameter  B is  unaffected  by  the 
modification. 


If  we  assume  a multiplicative  error  term  for  the  model  of 
Equation  (2.3.1),  the  model  of  observations  for  the  Arrhenius  model 
is  the  same  as  the  model  defined  in  the  previous  section  for  the 


power  rule  except  that  , g , x^  are  defined  as 


X.  = t 1 ’ (T± 1-T) 1 , i=1 k 


B'  = [B0,ei]  = [A*,B] 


xi  = (Ti1_T)  > 


i=l , . . . ,k  . 


(2.3.3) 


2 
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4 . Inference  at  Use  Stress  for  the 
Power  Rule  and  Arrhenius  Models 

The  inference  procedures  presented  in  this  section  apply,  in 
essence,  to  both  the  power  rule  and  tiie  Arrhenius  model. 

The  log  mean  life  at  use  stress,  log  0 , is  a linear  function 

° u 

of  the  unknown  parameters  6'  = > t.e., 

log  0 = 6.  + B.x 

n u 0 1 u 

where 

x * - log  V /V 
u u 

for  the  power  rule  with  V given  by  Equation  (2.2.3)  and 


(2.4.1) 


(2.4.2) 


ft  — 


x = T - T 
u u 


(2.4.3) 


for  the  Arrhenius  model  with  T given  by  Equation  (2.3.2).  In  Equation 


(2.4.2)  and  (2.4.3),  V and  Tu  are  use  stress,  respectively.  We 


obtain  the  best  linear  unbiased  estimate  of  log  0 , denoted  by 


log  0 , by  substituting  the  least  squares  estimates  in  Equation  (2.4.1) 


Thus , 


108  °u  = e0  + P’X 


1 A ’ 
1 11 


where  B^  and  Bj  are  the  least  squares  estimates  for  either  the  power 


rule  or  the  Arrhenius  model.  The  variance  of  log  ^ is  given  by 


Var(log  0 ) = Var(B^)  + x^Var(Bj) 


(2.4.4) 


= X (x'w_1x)  ]x’  , 

. u . u 


where  the  bottom  line  of  Equation  (2.4.4)  gives  the  variance  in  matrix 
notation  with  X = [1  , x ] . 


J 


T-353 


- 11  - 

The  exact  distributions  of  P^  and  P^  would  be  difficult  to 

obtain  since  they  are  linear  combinations  of  k independent  non- 
identically  distributed  log  gamma  random  variables.  The  log  gamma, 
however,  is  known  to  be  approxmately  normal  when  its  shape  parameter 
is  sufficiently  large.  Since  the  shape  parameter  of  z^  is  equal  to 

r^  , the  number  of  failures  observed  at  the  ith  stress  level,  we  may 

consider  the  estimates  to  be  linear  combinations  of  approximately  nor- 
mal random  variables  when  the  number  of  failures  observed  at  each  stress 
is  sufficiently  large.  Thus,  we  may  reasonably  assume  that  3q  and  3^ 

are  approximately  normal  if  the  number  of  failures  observed  at  each 
stress  level  is  moderately  large.  Also,  it  has  been  verified  that  Bq 

and  Bj  satisfy  Lyapunov's  condition  (see,  e.g.,  Parzen  [17],  page  432) 

for  convergence  to  normality  so  that  for  large  k 


~ ®1  ~ ^1 

^ N(0, 1)  and  - ^ N(0,1)  , 


Var(BQ) 


Var^) 


from  which  it  follows  that 


log  0 - log  0 log  0 - log  0 

■ -U - = U ■ U i N(0,1)  . 


Var(log  0 ) 


X (x'w-1x)  1x' 

~u  _ _ . _u 


Thus,  a large  sample  approximate  (l-a)%  confidence  interval  for  log  0 


is  given  by 


108  °u  ± Za/2VXu(X'“‘lxrl,<;  • 


(2.4.5) 


where  is  t^ie  a /2  percentage  point  of  the  standard  normal  dis- 

tribution. 

In  a typical  accelerated  life  test,  the  number  of  stress  levels 
k is  small  and  hence  the  degree  to  which  the  normal  approximates  the 
distribution  of  the  estimators  may  be  of  concern.  However,  even  for  a 
small  number  of  stress  levels,  the  estimators  are  approximately  normal 


T-  353 


- 12  - 

in  the  sense  that  their  skewness  and  kurtosis,  or  shape  factors,  are 
close  to  zero,  the  values  for  the  normal  distribution,  provided  a 
moderate  number  of  failures  is  observed  at  each  stress.  Also,  the  shape 
factors  approach  zero  as  the  number  of  failures  observed  and  the  number 
of  stress  levels  increases.  This  is  shown  in  Appendix  A. 

We  note  that  if  it  is  assumed  that  log  0 'v  N(log  0^  , 

X (X'W"1X)"1X’)  , it  follows  that  0 ^ A(log  0 , X (X 'W_1X) _1 X ’ ) 

.11  - -U  u u u ..u 

where  ^ A denotes  "approximately  log  normally  distributed."  Thus 
E(o"  ) ^ 0 exptX  (X'W_1X)_1X'/2]  , and  0 =0  exp[-X  ( X ' W~ 1 X) “ 1 X ’ / 2 ] 

U U .U  ^ . U U U - u .u 

is  an  approximately  unbiased  estimate  of  0 with 

u 

Var (0  ) = 02{exp[X  (X'W_1X)-1X' ] -1 } . 
u u u _ . u 

5 . Example  for  the  Power  Rule 

In  this  section,  we  illustrate  procedures  presented  in  previous 
sections  using  data  simulated  for  the  power  rule  by  Singpurwalla  [19].  An 
example  involving  the  Arrhenius  model  would,  of  course,  be  similar.  The 
data  was  generated  on  a computer  by  choosing  C = 1000  , P = 3 and 
-p 

using  0.  = CV  to  calculate  five  values  for  . . The  results  of 

it  i 

the  simulation  are  shown  in  Table  2.1  with  the  n.  and  r.  used  for 

l l 

each  and  0^  determined  by  Equation  (2.1).  Also  shown  are  the 

approximately  unbiased  estimate  of  0.  , and  the  8 and  V(8)  deter- 
mined by  Equations  (2.2,5)  and  (2.2.6),  respectively. 

The  maximum  likelihood  estimate  of  P , 3.09,  is  comparable  to 
the  least  squares  estimate  of  P , 6^  = 3.05  . The  variance  of  the 

least  squares  estimate  is  2.94%  greater  than  the  Cramer-Rao  lower  bound 
on  the  variance  of  the  maximum  likelihood  estimate.  The  estimates  of 
C*  are  not  comparable  because  the  value  of  C*  in  the  least  squares 
case  is  different  from  the  value  of  C*  in  the  maximum  likelihood  case. 
That  is,  the  modification  of  the  power  rule  required  to  obtain  uncor- 
related least  squares  estimates  results  in  a different  value  of  C* 


k 


TABLE  2.1 

RESULTS  OF  A MONTE  CARLO  ACCELERATED  LIFE  TEST 


n, 

r . 

0. 

0J 

i 

l 

l 

i 

l 

i 

10 

30 

15 

1.000 

1.308 

0.998 

20 

30 

15 

0.125 

0.078 

0.123 

30 

30 

20 

0.037 

0.030 

0.036 

40 

30 

25 

0.016 

0.017 

0.015 

50  ' 

30 

25 

0.008 

0.008 

0.008 

Note : 

With 

the  ' 

exception  of 

0^  , this 

data  ii 

from  Singpurwalla  [19]. 


Least  Squares  Estimates  for  the  Power  Rule  Based  on  the  Above  Data: 


8'  = (-3.2526  , 3.0464)  . 


Variance-covariance  Matrix  of  the  Least  Squares  Estimates: 


'0.01025  0 

Var( 8)  = 

0 0.03538 
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from  that  which  results  when  the  power  rule  is  modified  to  yield 
asymptotically  uncorrelated  maximum  likelihood  estimates.  Also,  the 
maximum  likelihood  procedure  estimates  C*  while  the  least  squares 
estimates  log  C*  . For  = 7 , the  least  squares  estimate  of  log  0 

is  log  0 = 1.1105  and  the  approximately  unbiased  estimate  of  0 
u u 

is  0 = 2.9126  . These  compare  to  the  actual  values  of  log  0 = 1.0700 

u u 

and  0 = 2.9155  . 

u 

The  skewness  and  kurtosis  of  the  least  squares  estimates 

and  3.  obtained  from  the  data  in  Table  2.1  are  y ' = -0.1012  , 

1 ,30 

y„  " = 0.0205  and  y.  ' = -0.1552  , y " = 0.0549  . These  values 

->0^  1 ,Pj  2,3j 

are  reasonably  close  to  zero  so  that  we  may  assume  and  6^  are 

approximately  normal  and  it  follows  that  log  0 is  approximately 

normal.  Thus,  Expression  (2.4.5)  may  be  used  to  construct  approximate 
(l-a)%  confidence  intervals  for  log  0 . With  Var(log  0^)  = 0.0828  , 

an  approximate  95%  confidence  interval  for  log  0 is  obtained  as 

1.1105  + 1.96/0.0828  = [0.5465,  1.6745]  . 

6 . Eyr ing  Model 

In  this  section  we  consider  the  Eyring  model  for  an  exponential 
scale  parameter.  In  a recent  paper,  Singpurwalla,  Castellino  and 
Goldschen  [22],  obtained  least  squares  estimates  for  the  parameters  of 
this  model.  In  this  section,  we  present  a modification  of  the  Eyring 
model  that  results  in  orthogonal  least  squares  estimators  for  its  param- 
eters. Since  the  regression  model  defined  in  [22]  results  in  an  ill- 
conditioned  covariance  matrix,  an  orthogonal  design  for  the  estimation 
of  the  Eyring  parameters  is  desirable. 

An  exponential  distribution  was  assumed  in  [22]  and  a generalized 

form  of  the  Eyring  model  was  used  to  re-parameter ize  the  scale  parameter 

A.  , in  terms  of  a thermal  stress  T.  and  non-thermal  stress  V.  . 
i t l 

The  Eyring  model  assumes  that  the  time  rate  degradation  of  A.  may  be 


I 
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expressed  in  terms  of  the  stress  pair  (T^,V^)  as 


A,  = AT.  exp  [-B/KT.  + CV . + DV./KT.] 
i i x l li 


where  K is  Boltzmann's  constant  and  A,  B,  C,  D are  unknown  parameters. 

Estimates  of  the  unknown  parameters  may  be  obtained  from  censored 
sample  life  tests  conducted  at  k > A different  stress  environments 
{T.,V.}  , i = l,...,k.  The  ith  life  test  is  conducted  as  described 

in  the  first  section  of  this  chapter.  An  estimate  of  the  hazard  rate 
A.  may  then  be  obtained  by  taking  the  reciprocal  of  the  estimate 

0 , i.e.  A.  = 0,^  . Since  0.  has  a gamma  distribution.  A,  has  an 
l i i i l 

inverted  gamma  distribution  with  scale  parameter  and  shape 

parameter  r^  . 

We  may  now  formulate  a linear  model  for  the  Eyring  parameters 

"-1 

in  terms  of  log  X.  . Since  A.  = 0 , log  A.  = -log  0.  and  since 

°i  il  l i 

log  A^  is  merely  the  negative  of  log  0^  , the  two  variables  have 
similar  properties.  Thus,  log  A^  has  a log  gamma  distribution  with 
mean  and  variance 


E(log  A1)  = log  Aa  + log  - <Kri) 


(2.6.1) 


Var (log  AJ  = <|r  (r^  . 


(2.6.2) 


Orthogonal  least  squares  estimates  for  the  Eyring  parameters  may  be 
obtained  by  modifying  the  original  form  of  the  Eyring  model.  The  basic 
character  of  the  modified  model  is  the  same  as  the  original  and  is  given 


1 ‘ A*Tl  e4-B(icT7  -T)  + C<VV)  + d(k^  - VT)J 


where 


- 16  - 
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k w.  (KT.) 

y — J »L_ 


(2.6.4) 


k w.1  V. 


V = Z 


j=i  ; -i 
J E w . 

3 = 1 3 


(2.6.5) 


k w.  (V./KT.)  1 

VT  = £ - 

3=1  J -1 
E w . 

3=1  3 


(2.6.6) 


w.  = <P  (r  .)  , j = 1, . . . ,k 


(2.6.7) 


In  effect,  the  modification  is  a change  in  the  proportionality  constant 
indicated  by  the  notation  A*  . The  values  of  B,  C and  D are  unaffected 
by  the  modification. 

If  we  assume  a multiplicative  error  term  for  the  model  given  by 
Equation  (2.6.3)  a linear  model  of  observation  for  the  Eyring  model  can 
be  written  as 

Z = X B + e 


where 


X = • 


kit. 


T-  3 5 3 


X. 
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/ 1 -\ 

'(kvt) 

, (V.-V) 

B1,B2,B3,B 

] - [log 

[log  ex,... 

•los  ek] 

[zl zk] 

i = 1 k 


, -B,C,Dj 


zi  = log  - log  rt  + iKr±)  - log  T.  , i = l,...,k  . 

From  Equations  (2.6.1),  (2.6.2)  and  the  randomization  of  the  life  test 

E(e)  = 0 

E(e'e)  = W = diag[w1 , . . . ,wk]  . 

The  w.  are  given  by  Equation  (2.6.7)  and  W is  a k x k diagonal 

matrix  which  is  the  variance-covariance  matrix  of  e . 

Equations  for  the  least  squares  or  best  linear  unbiased  estimate 
of  8,3  and  Var(3)  , are  given  in  matrix  form  by  Equations  (2.2.4) 
and  (2.2.5).  For  the  modified  form  of  the  Eyring  model  these  equations 
reduce  to 


B " ( Bl, B2. ^3, Ba) 

[Ew.^z.  Ex.-w.^z. 

li,  i2  l l 

_ -1  r 2 -1 
Ew . Ex . _w. 

1 i2  i 


Ex 

Ex 


-1  _ -1 
■ z . Ex . . w . z . I 
i3  i i , i4  l i| 

2-1  _ 2 -1 
x . _w.  Ex . . w . J 

l3  1 14  1 


(2.6.8) 


and 


Var(8)  = diagpar(81),  Var(82),  Var(B3),  Var(B4)J 

= diag^Ew^)  1,  (Ex^2w11)  \ (Ex^w^)  1,  (Ex^w^j 


where 
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where  the  second  line  gives  the  variance  in  matrix  notation  with 

X = [1  x0  X-  X ] . 

2u  3u  4u 

The  situation  in  this  case  is  similar  to  that  considered  previously 
for  the  power  rule  and  Arrhenius  models.  That  is,  the  exact  distribu- 
tion of  the  estimates  B^  , i = 1,...,4  would  be  difficult  to  obtain 

since  they  are  linear  combinations  of  non-identical  log  gamma  random 
variables.  However,  the  log  gamma  may  be  considered  to  be  approximately 
normal  if  the  gamma  shape  parameter  is  sufficiently  large.  Thus,  if  the 
number  of  failures  observed  at  each  stress  level  is  sufficiently  large, 
we  may  consider  the  estimates  to  be  linear  combinations  of  approximately 
normal  random  variables.  Also,  it  has  been  verified  that 
B^  , i = 1,...,4  satisfy  Lyapunov's  condition  for  convergence  to 

normality.  Thus,  for  large  k 


B.  - 6. 

— j— * N(0,1)  , i = 1,2, 3, 4 
fvar(B.) 


and  it  follows  that 


log  Xu  - log  Au 


X 8 - X 3 

■.LK  .IK 


VVar ( log  A ) J X (X,W_1X)~1X' 
u T u ^ „u 


-v  N(0, 1) 


Thus,  a large  sample,  approximate  (l-a)%  confidence  interval  for 
log  A is  given  by 


log  A + Za/2  WVar(log  A ) 
u — » u 


= X 6 + log  T + Za/2  V*  (X'W_1X)  V 
-U.  u — * .u  - . .u 


where  Za/2  is  the  a/2  percent  point  of  the  standard  normal  distribu- 


tion . 
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In  a typical  accelerated  life  test,  the  number  of  stress  levels 
k is  small  and  hence  the  degree  to  which  the  normal  approximates  the 
distribution  of  the  estimates  may  be  of  concern.  We  may,  however, 
obtain  estimates  that  are  approximately  normal  even  for  a small  number 
of  stress  levels.  If  the  number  of  failures  observed  at  each  stress  is 
moderately  large,  the  estimates  are  approximately  normal  in  the  sense 
that  their  skewness  and  kurtosis,  or  shape  factors,  are  close  to  zero, 
the  values  for  the  normal  distribution.  Also,  the  shape  factors  of  the 
estimates  approach  zero  as  the  number  of  failures  observed  and  the  num- 
ber of  stress  levels  increase.  This  is  shown  in  Appendix  B. 

We  note  that  when  the  best  linear  unbiased  estimate  of  log  A 

is  approximately  normally  distributed,  i.e., 

— • -1-1 
log  A 'v  N(log  X , X (X’W  X)  X')  , it  follows  that 
u u u - n 

— • -1-1 

A 'v  A(log  A , X (X’W  X)  X1)  where  A"  denotes  "approximately  log 
u u n u 

normally  distributed."  Thus,  E(A  ) ^ A exp[X  (X’W  ^X)  ^ X ' / 2 ] and  an 

u u _ u _ _ u 

approximately  unbiased  estimate  of  A^  is  given  by 

A = A exp [ -X  (X'W_1X)"1X'/2]  with  Var(A  ) = A (exp[X  (X'W-1X)-1X’ ] -1 } 
u u u u u u u .......  u 

8.  Example  for  the  Eyring  Model 

Some  accelerated  life  test  data  for  test  items  which  were  sub- 
jected to  both  temperature  and  voltage  stress  is  given  in  Table  2.2. 

This  data  was  published  in  [22]  and  is  used  here  to  illustrate  estima- 
tion procedures  for  the  Eyring  model.  Orthogonal  least  squares  estimates 
for  the  Eyring  model  parameters  based  on  this  data  are  shown  in  Table  2.3 
along  with  the  variances  and  shape  factors.  The  equations  for  the  shape 
factors,  or  skewness  and  kurtosis,  of  the  estimates  for  the  Eyring 
parameters  are  derived  in  Appendix  B.  The  values  of  the  shape  factors 
for  the  estimates,  shown  in  Table  2.3,  are  fairly  close  to  zero  so  that 
an  assumption  of  normality  appears  to  be  reasonable. 

Use  stress  values  for  the  items  used  to  obtain  the  life  test 
data  of  Table  2.2  were  not  available.  However,  if  we  assume  that  the 
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estimates  are  approximately  normal  then  log  \ or  any  other  linear 
combination  of  the  estimates  would  also  be  approximately  normal. 

The  accelerated  life  test  data  shown  in  Table  2.2  illustrate  a 
point  in  regard  to  measuring  the  approximate  normality  of  the  Eyring 
parameter  estimates  by  their  shape  factors.  That  is,  since  the  shape 
factors  of  the  estimates  are  weighted  averages  of  the  shape  factors 
of  the  individual  observations,  the  effect  of  obtaining  a small  number 
of  failures  at  one  or  more  stress  levels  may  be  offset  by  observing 
a larger  number  of  failures  at  other  stress  levels.  Thus,  the  shape 
factors  of  the  estimates  obtained  from  the  data  in  Table  2.2  approxi- 
mate the  values  for  the  normal  distribution  although  only  one  failure 
was  observed  at  two  of  the  stress  levels. 
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TABLE  2.2 

TANTALUM  CAPACITOR  FAILURE  TIMES  UNDER 
VOLTAGE  AND  TEMPERATURE  STRESSES 


Test 
Number 
[Index  i) 

Number 
on  Test 

n . 

i 

Vo  I tage 
V. 

l 

Tempera- 
ture °c 
T. 

l 

Actual  Times 
to  Failure 
(hr) 

Number  of 
Failures 
r . 

l 

1 

1,000 

35.0V 

85  °C 

20,90,700, 

37000 

4 

2 

200 

40.6V 

85  °C 

20,3600,9500, 

27000 

4 

3 

50 

46.5V 

85  °C 

800,2800 

2 

4 

50 

51.5V 

85  °C 

500,800,2400, 
10700,  10700 

5 

5 

500 

46.5V 

45  °C 

100,1200,7500, 

20000,26000, 

27300 

6 

6 

175 

46.5V 

5 °C 

1000 

1 

7 

175 

62.5V 

5 °C 

25,40,165,500, 

620,720,820, 

910,980,1270, 

1600,2270, 

2370,4590, 

4690,4880. 

7560,8730, 

12500 

19 

8 

50 

57.0V 

45  °C 

8900 

i 

Data  from  Singpurwal la , Castellino  and  Goldschen  [22], 

TABLE  2.3 


ORTHOGONAL  LEAST  SQUARES  ESTIMATES  FOR 
THE  EYRING  MODEL  PARAMETERS  BASED 
ON  THE  TEST  DATA  OF  TABLE  2.2 


i 

fl. 

1 

. 

Var (6 . ) 
i 

v 

Y1  ,Bt 

Y2,P. 

1 

i 

-18.4933 

0.0261 

-0.1604 

0.0512 

2 

0.4054xl0-12 

0. 3729xlO_26 

-0.1789 

0.0688 

3 

0.1543 

0.2709xl0-3 

-0.2032 

0.0926 

4 

0.3870xl0"U 

-30 

0.2062x10 

-0.1831 

0.0715 
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CHAPTER  III 


A BAYESIAN  ANALYSIS  OF  THE  POWER  RULE  MODEL 


1.  Introduction 

In  this  chapter  we  consider  a Bayesian  approach  for  the  analysis 
of  accelerated  life  test  data.  Here  the  power  rule  is  used  to  re- 
parameterize the  scale  parameter  of  an  exponential  and  a Weibull  dis- 
tribution. The  approach  we  take  is  the  following:  Re-parameterize 

the  distribution  of  the  estimator  of  the  scale  parameter  in  terms  of 
the  power  rule  parameters.  Specify  an  appropriate  joint  prior  distri- 
bution and  compute  the  posterior  distribution  using  data  from  accel- 
erated life  tests.  The  joint  posterior  distribution  then  provides  a 
solution  to  the  problem  of  inference  concerning  the  unknown  parameters. 


2.  The  Exponential  Case 

The  power  rule  model  for  use  with  the  scale  parameter  of  an  ex- 
ponential distribution  was  described  in  Chapter  II.  Accelerated  life 
tests,  conducted  as  described  in  Chapter  II,  provide  estimates  of  0^ 

to  be  used  for  a Bayesian  analysis  of  the  power  rule.  The  joint  like- 
lihood of  0 = (0  ,...,0^)  given  the  estimates  0 = (§,..., §^)  ob- 
tained from  tests  conducted  at  k distinct  accelerated  stress  levels 
is  based  on  the  distribution  of  0^  , i=l,2,...,k  , and  is  given  by 


L(0/0)  = 


k 

n 

Li=l 


r r -1 
(ri/01)  1(§i) 


ro^) 


A likelihood  in  terms  of  C and  P may  be  obtained  by  substituting 

p 

0.  = C/V^  , i=l,...,k  in  the  above  expression.  A more  convenient 
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parameterization  is  obtained  by  slightly  modifying  the  power  rule  to 


where 


H.  = C/(V./V)P  , 
I 1 


(3.2.1) 


This  modification  results  in  a different  value  for  C (in  effect,  a 
change  in  the  proportionality  constant),  while  leaving  the  value  of  P 
unaffected.  The  likelihood  is  then 


L(C,P/0)  = 


k 

n 

i=l 


„iri 


r.(V1/V) 


(0t) 


V1 


r(r.) 


exp 


k r .0 . (V. /V)P 

' I ^ 

i=l  L 


(3.2.2) 

In  this  case  the  modified  power  rule  gives  the  likelihood  a more  circu- 

p 

lar  shape  than  would  be  obtained  using  0^  = C/V^  and  also  has  the 


effect  of  making  the  likelihood  easier  to  manage  computationally. 


The  usual  Bayesian  approach  requires  that  a prior  distribution 
n(C,P)  be  specified  in  order  to  express  a state  of  knowledge  regarding 

A A A 

C and  P before  data  is  obtained.  Given  the  data  0 = (9  6^) 

and  a prior  distribution  I1(C,P)  , a posterior  distribution  for  (C,P) 
may  be  obtained  from  the  relationship 

p(C,P/0)  « L(C, P/0)  n(C,P)  . 

The  viewpoint  held  by  many  Bayesians  (e.g..  Box  and  Tiao  [3],  p.  20) 
maintains  that  a statistical  inference  problem  is  regarded  as  solved 
when. an  appropriate  probability  statement  can  be  made  regarding  the 
state  of  nature  (i.e.,  the  unknown  parameters)  in  question.  Thus,  the 
posterior  distribution  p(C,P/§)  provides  the  solution  to  the  problem 
of  what  may  be  inferred  about  the  parameters  C and  P from  the  data 
0 and  a relevant  state  of  prior  knowledge  represented  by  I1(C,P)  . 
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It  would  be  reasonable  to  assign  a continuous  prior  distribution 
for  C and  P , since  C and  P are  defined  on  a continuous  param- 
eter space.  This  presents  difficulties  in  the  sense  that  the  integrals 
that  must  be  evaluated  in  order  to  normalize  the  posterior  distribution 
do  not  possess  convenient  closed  form  solutions.  In  such  cases,  it  is 
usual  to  resort  to  numerical  integration.  As  an  alternative  to  numeri- 
cal integration  we  consider  a simpler  numerical  procedure.  'Phis  proce- 
dure involves  a discretization  of  the  prior  distribution. 

For  our  problem,  this  procedure  assumes  that  prior  knowledge 
regarding  C and  P may  be  expressed  in  terms  of  a discrete  bivariate 
probability  mass  function.  The  process  of  obtaining  a discrete  prior 
for  (C, P)  may  be  considered  consistent  with  usual  Bayesian  procedures 
in  the  sense  that  the  prior  may,  in  effect,  be  specified  before  observ- 
ing the  data.  For  instance,  a continuous  distribution  could  be  selected 
as  representing  prior  knowledge  regarding  (C,P)  and  then  appropriately 
discretized.  For  computational  efficiency  it  is  most  economical  to 
discretize  the  prior  only  over  certain  regions  of  the  parameter  space 
which  are  dictated  by  the  available  data.  That  is,  after  observing  the 
data  a region  R of  the  parameter  space  is  identified  for  which  the 
likelihood  is  appreciably  non-zero  and  the  prior  discretized  over  R . 
The  general  form  of  R is 

R = { (C,  P) : a < C < B,  Y<.P<6}, 

where  a,  6.  Y.  6 are  boundary  points  such  that  outside  R the  like- 
lihood is  negligible.  The  subset  of  the  parameter  space  outside  R is 

therefore  not  relevant  to  the  realized  likelihood  and  thus  little  or 
nothing  is  lost  by  excluding  it  from  consideration. 

After  R has  been  identified,  a discrete  prior  may  be  obtained 
by  dividing  R into  rectangular  subsets  and  assigning  probability  mass 
to  each  subset.  This  requires  that  intervals  for  C and  P , say 

and  » respectively,  be  chosen  so  that  the  subsets  defined  by 

R = { (C,P) : a+rnA,  < C < a+(nH-l)A,  , Y+nA„  < P < Y+(n+l)A.}  , 

mn  1 1 Z Z 


m=0, 1, . . . , M , n=0, 1 , . . . ,N  , 
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provide  a suitable  partition  of  R . The  probability  mass  assigned  to 

R , denoted  bv  W , is  then  associated  with  the  mid-point  of  R , 
mn  mn  mn 

(C  = a+(m+7r)A,,  P = Y+(n+^-) A_ ) . The  discrete  prior  distribution  is 
m l 1 n l L 

then 


II(C  ,P  ) = W 
m n mn 


m=0,l,...,M  , n=0,l,...,N  , 


where  fl ( C ,P  ) may  be  interpreted  as  the  approximate  prior  probability 
m n 

that  C = C and  P = P . The  posterior  is  then 
m n 


P (C  ,P  /6)  = 
m n ~ 


l(c  ,p  /e>n(c  ,p  ) 


m n 


m=0,l,...,M  , n=0,l,...,N  , 


where  the  likelihood  L(C  ,P  /0)  is  given  by  Equation  (3.2.2),  and 

m n ~ 


M N 

F-  l l L(C„V?> 

m=0  n=0 


It  is  difficult  to  specify  what  constitutes  a suitable  partition 
of  R . Generally,  the  interval  lengths  A^  and  A2  should  be  of  a 

size  to  adequately  represent  the  behavior  of  the  likelihood  over  R . 

If  the  intervals  are  too  large,  much  of  the  variation  in  the  posterior 
will  be  omitted,  and  if  they  are  too  small,  the  computations  will  be- 
come unwieldy. 


After  R and  A^,A2  have  been  determined,  obtaining  a discrete 

prior  distribution  based  on  a continuous  distribution  is  reasonably 
straightforward.  Suppose  f(x,y)  is  the  p.d.f.  of  an  absolutely  con- 
tinuous distribution  which  is  an  appropriate  representation  of  prior 
knowledge  regarding  (C,P)  . The  probability  masses  for  a discrete 
prior  distribution  may  then  be  obtained  by  evaluating  the  probabilities 


Pr[a+mAj  <_  X <_  a+(nri-l)A^,  < Y < y+(n+l)A2]  = 


a+(nri-l)A 

y+(n+l)A2 

f 

/ f(x,y)dydx  = W 

a+mAj 

y-^a2 

m=0,l,...,M  , n=0,l,...,N  . 


k 


An  alternative  method  for  constructing  a discrete  prior  for 
(C,P)  is  to  subjectively  assign  probability  masses  to  subsets  of  the 
parameter  space.  After  the  region  R has  been  identified,  we  may  as- 
sign a subjectively  chosen  probability  mass  W to  each  of  the  subsets 

R . For  example,  a discrete  locally  uniform  prior  over  the  region  R 
mn 

is  obtained  by  assigning  a constant  probability  to  each  of  the  subsets 

R . The  procedure  is  quite  general  in  the  sense  that  any  values  for 
mn 

the  W may  be  assigned, 
mn 

• 2.1  Example 

The  results  of  a simulated  accelerated  life  test  were  given  in 
Table  2.1.  In  this  section  we  use  this  data  to  illustrate  the  Bayesian 
methods  presented  in  this  chapter. 


We  make  the  assumption  here  that  prior  knowledge  regarding  (C,P) 
may  be  expressed  in  terms  of  a joint  prior  distribution  which  assigns  a 
positive  probability  for  non-negative  values  of  C and  P . A joint 
distribution  which  satisfies  this  requirement  is  the  absolutely  contin- 
uous bivariate  exponential  (ACBVE)  distribution  given  by  Block  and  Basu 
[2].  The  ACBVE  is  defined  for  0<x<°°,  0<y<°°.  In  order  to 
have  more  flexibility  in  specifying  the  prior,  especially  in  connection 
with  the  location  of  the  observed  likelihood,  we  introduce  location  pa- 
rameters y^  and  ^2  • The  density  for  the  shifted  form  of  the  ACBVE 
is  then 


^1^^2+^12^e 


•Ajfx-Mj)  - (X2+X12) (y-h2) 


-yvV  -yw 

2A  - (A2+A12)e  1 (A^A^je  Z 1 Z 


f(x,y)  = 


y ^i+^y 6 


for  y^  < x < y 


-^i+*i2^x_y  ~ yy-y 


-A  (y.-y  ) -A  (y  -y  ) 

2A  - (A2+A12)e  1 (Ayye 


for  y2  < y < x 


..... 
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A = Ai  + X2  + A12 


,A2’X12 


> 0 


VU2 


> 0 


By  introducing  the  location  parameters  we  also  avoid  assigning  high 
prior  probability  near  the  origin. 


Three  contours  of  posterior  distributions,  together  with  the 

A A 

mode  (C,P)  , obtained  using  different  prior  distributions  are  shown 
in  Figure  3.1(i)  - (vii).  The  levels  of  these  contours  are  defined  by 
the  following: 


A:  p (C, P/0)  >0.5  p(C,P/§) 
B:  p(C,P/6)  > 0. 25p(C, P/0) 
C:  p(C,P/6)  > 0.05p(C,P/6) 


That  is,  the  probability  associated  with  values  of  (C,P)  within  the 
A contour  is  at  least  0.5  as  large  as  the  probability  associated  with 
the  mode;  within  B,  at  least  0.75  as  large  as  the  mode;  and  within  C, 
at  least  0.95  as  large  as  the  mode.  Using  the  bivariate  normal  approx- 
imation, these  contours  correspond  roughly  to  the  boundaries  of  the 
50,  75,  and  95%  highest  posterior  density  (H.P.D.)  regions,  respectively 
(see  Box  and  Tiao  [3],  pp.  122-125). 

The  posterior  contours  for  a uniform  prior  are  shown  in  Figure 
3.1(i).  The  posterior  contours  shown  in  Figure  3.1(ii)  - (vii)  were 
obtained  using  ACBVE  priors  with  parameter  values  given  in  Table  3.1. 

In  each  case,  a = 0.0285,  y = 2.225,  Aj  = 0.001,  = 0.05,  M = 29, 

and  N = A3  . The  plots  in  Figure  3 . 1 ( i ) - (vii)  are  similar  so  that 
for  the  cases  considered,  the  priors  have  little  effect  on  the  posterior 
contours.  Also,  for  each  case,  the  mode  (C,P)  is  (0.039,  3.05)  . 

The  similarity  of  these  results  may  be  accounted  for  by  the  dominance 
of  the  realized  likelihood  in  a small  region  of  the  parameter  space 
and  the  lack  of  variation  of  the  ACBVE  priors  in  this  region.  We  em- 
phasize here  that  these  results  apply  only  for  the  data  analyzed  and 
that  another  set  of  data  would  yield  different  results. 


4 «' 


TABLE  3.  I 


PARAMETER  VALUES  FOR  THE  BIVARIATE  EXPONENTIAL 
PRIORS  USED  TO  OBTAIN  THE  POSTERIOR  CONTOURS 
SHOWN  IN  FIGURE  3.1(i)  - (vii) 


A1 

| 

; 

A12 

U1 

M2 

ii) 

25.0 

0.20 

0. 18 

0.0285 

2.2: 

ii) 

15.0 

0.25 

0.20 

0.0285 

2.2: 

iv) 

0.  5 

0.08 

0.  10 

0.0285 

2.2: 

(v) 

0.5 

0.08 

0. 10 

0.0 

2.2: 

vi) 

0.5 

0.10 

0.0285 

0.0 

.5 
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For  comparison,  the  exact  50,  75,  and  95%  contours  and  the  mode 
of  the  bivariate  normal  distribution  corresponding  to  this  problem  are 
shown  in  Figure  3.1(viii).  Following  a theorem  of  Lindley  ([10],  p. 
130),  posterior  distributions  are  asymptotically  approximately  multi- 
variate normal  with  means  equal  to  the  maximum  likelihood  estimates 
and  variance-covariance  matrix  equal  to  the  information  matrix  evalu- 
ated at  the  maximum  likelihood  estimates.  Strictly  stated,  Lindley' s 
theorem  does  not  apply  in  this  case  since  we  have  a priori  restricted 
(C,P)  to  positive  values.  However,  the  bivariate  normal  contours 
based  on  the  maximum  likelihood  estimates  for  (C,P)  are  such  that 
negative  values  have  negligible  probability.  The  contours  of  the  pos- 
terior distributions  are  roughly  similar  in  shape  and  location  to  the 
bivariate  normal  contours.  The  asymptotic  normality  of  the  posterior 
depends  on  the  sample  size  which  in  this  case  corresponds  to  the  number 
of  stress  levels  and  the  number  of  failures  per  stress  level.  Thus,  as 
the  number  of  stress  levels  and  the  number  of  failures  per  stress  level 
increase,  we  would  expect  the  agreement  of  the  posterior  contours  with 
the  bivariate  normal  contours  to  be  enhanced. 


2.2  Prediction  at  Use  Stress 

Since  we  assume  that  the  power  rule  relationship  holds  at  use 

stress  conditions,  a particular  predicted  value  for  0 is  obtained  as 

u 


-P 


0 , , - C (V  /V)  " , 

u(m,n)  m u 

where  V is  use  stress  and  V is  given  by  Equation  (3.2.1).  Given 


the  discrete  posterior  distribution  for  (C,P)  , the  probability  asso- 
ciated with  0 . . is 

u(m,n) 


p 1 0 , , = c (v  /v)  n / e I = P(c  ,p  /§)  . 

|_  u(m,n)  mu  ' ~j  mn~ 

Thus,  a predictive  distribution  for  0^  may  be  obtained  based  on  the 

posterior  distribution  of  (C,P)  . It  also  follows  that  an  approximate 
95%  highest  posterior  density  (H.P.D.)  region  for  0^  may  be  determined 


■Hi 
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based  on  the  95%  H.P.D.  region  of  p(C  ,P  /0)  . That  is,  the  values 

m n 

of  0 , v for  (C  ,P  ) such  that  p(C  ,P  /0)  > 0.05  p(C,P/0)  con- 
u(m,n)  mn  m n ~ — 

stitute  the  set  of  0 values  within  the  approximate  95%  H.P.D.  region 
for  0 

u 

By  grouping  the  predicted  values  for  0^  within  the  approximate 

95%  H.P.D.  region  and  summing  the  probabilities  associated  with  each 
group,  we  may  obtain  a distribution  for  0 . In  order  to  arrange  the 

predicted  values  for  0 in  groups,  we  define 

6u(j)  = °u(0)  + jA/2  ’ j=l»2,...,J  , 

where  0 ...  is  the  lower  limit  for  values  of  0 within  the  approxi- 
u(0)  u 

mate  95%  H.P.D.  region,  A is  an  arbitrarily  selected  interval  size, 
and  J is  the  number  of  intervals.  We  note  that  J and  A are  se- 
lected so  that  the  largest  value  predicted  for  0 is  between 

u 

0 ...  + A/2  and  0 , T.  - A/2  . The  probability  associated  with 
u(J)  ll(J) 

0 , T.  is  obtained  as 
u(  J) 

P (0  ,.v/0)  = l P(0  f JQ)  , 

u(j)  ~ (m,n)  u(m,n) 


where 


(m,n) . = { (m,n) 


0 . . < 0 ...  + A/2  and  0 . . > 0 , . . - A/2}  , 

u(m,n)  — u(j)  u(m,n)  — u(j) 


j— 1,...,J  . 


Some  distributions  of  predicted  values  for  0 at  V = 7 are 

u u 

shown  in  Figure  3.2.  These  were  obtained  from  the  posterior  distribu- 
tions with  contours  shown  in  Figure  3 . 1 ( i ) - (iv).  The  labels  (i)  - 
(iv)  in  Figure  3.2  correspond  to  the  cases  shown  in  Figure  3.1(i)  - (iv). 
As  would  be  expectec}  from  the  similarity  of  the  posterior  contours  for 

(C,P)  , the  prediction  distributions  for  0 are  quite  similar. 

u 
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The  approximate  93%  probability  limits  for  0 are  obtained  as 

u 

the  largest  and  smallest  0 . . such  that  p(9  , . /§)  > 

u(m,n;  u(m,n)  ~ — 

0.05  p(C,P/§)  . The  approximate  95%  probability  limits  for  6 based 
~ u 

on  the  posterior  distributions  with  contours  shown  in  Figure  3.2  are 

shown  in  Table  3.2.  The  values  of  these  probability  limits,  for  the 

different  cases,  are  the  same  or  very  close. 

3.  The  Weibull  Case 

In  this  section  we  consider  the  use  of  the  power  rule  with  a two- 
parameter  Weibull  distribution.  The  cumulative  distribution  of  the  Wei- 
bull may  be  written  as 

- (t/0  )b 

F(t;  6i,b)  = 1 - e 1 , 

where  the  scale  parameter  0 is  assumed  to  be  a functi >n  of  the  ith 

stress  level.  The  Weibull  is  useful  as  a failure  distribution  due  to 
its  well  known  flexibility  in  reflecting  the  failure  characteristics  of 
components  and  devices. 

We  consider  the  power  rule  as  a re-parameterization  for  the  Wei- 
bull scale  parameter  in  the  following  form: 

= 1/CV*  , 

where  describes  the  stress  and  C and  P are  unknown  parameters. 

In  this  section  we  use  a Bayesian  approach  for  inference  concerning  the 
power  rule  parameters. 

When  the  shape  parameter  b is  known,  the  Weibull  c.d.f.  may  be 

written 

-x/n. 

F(x;  n1)  = 1 - e , 

b b 

where  x = t and  = 0^  . Thus,  x has  an  exponential  distribution 
with  scale  parameter  • An  estimate  of  may  be  obtained  from  life 
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tests  conducted  as  described  in  Section  2 of  this  chapter.  If  the  ith 
test  is  terminated  after  r.  (_<  n^)  failures  occur,  a minimum  variance, 

unbiased  efficient  and  sufficient  estimator  of  r).  is  (Epstein  and 

Sobel  [6]) 


b ^ . . b 

r t . . + (n.-r . )t . 
i il  it  lr. 

~ v i 

ni  = 7. 

j=i  i 


where  t , j=l, . . . 
stress  level.  The  d 
scale  parameter  (r^ 


,r.  are  the  observed  failure 

l 

istribution  of  ru  is  a gamma 
/n .)  and  shape  parameter  r. 


times  at  the  ith 
distribution  with 


The  approach  for  obtaining  posterior  distributions  for  the  power 
rule  parameters  in  this  case  parallels  that  of  the  previous  sections 
for  the  exponential  distribution.  In  constructing  a likelihood  in 
terms  of  C and  P we  modify  the  power  rule  to 

eA  = 1/C(V./V)P  , 


where  V is  given  by  Equation  (3.2.1).  In  this  case,  the  likelihood 
is  conditional  on  the  data  n = (n^ , . . . ,fj^)  and  the  value  of  b . That 
is. 


L(C,P/n,b)  = 


k 

n 

i=i 


t. 


cb(v/v')Pb 


3 i r.-l 
] ' 1 


r(r  ) 


• exp 


l r n.cb(v  /V)pb 
i=i 


The  use  of  the  modified  power  rule  in  this  case  results  in  the  likeli- 
hood being  computationally  easier  to  manage. 


A posterior  distribution  for  C and  P may  be  obtained  using 
the  relationship  p(C,P/n,b)  <*  L(C, P/fj,b)Il(C,P)  , where  FI (C , P)  is  an 
appropriate  prior  distribution.  Using  the  discrete  procedure,  the  pos- 
terior is 


p(C 


,P  /U,b)  = 
n ~ 


L(c  ,p  /n,b>n(c  , 

m n ~ m 


P ) 

n 


m=0,l,...,M  , n=0,l,...,N 
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where  II(C  ,P  ) is  an  appropriately  discretized  prior  distribution 
m n 

M N 

and  F = T y L(C  ,P  /fj,b)II(C  ,P  ) . Procedures  for  obtaining 
„ ,,  m n ~ m n 

m=0  n=0 

FI ( C ,P  ) and  its  interpretation  given  in  previous  sections  dealing 
in  n 

with  the  exponential  case  also  apply  in  the  Weibull  case  under  consid- 
eration here. 

3 . 1 Examp  1 e 

Some  accelerated  life  test  results  published  by  Nelson  [15]  are 
shown  in  Table  3.3.  The  data  consists  of  breakdown  times  of  an  insu- 
lating fluid  under  different  voltage  stresses.  We  assume  that  the  fail- 
ure times  have  a Weibull  distribution  with  scale  parameter  0^  re- 
parameterized in  terms  of  the  stress  by  the  power  rule,  and  shape  pa- 
rameter b independent  of  the  stress. 

We  assume  here,  as  with  the  exponential  case,  that  prior  knowl- 
edge regarding  (C,P)  may  be  expressed  in  terms  of  a joint  prior  which 
assigns  positive  probability  for  non-negative  values  of  C and  P . 

Thus,  the  absolutely  continuous  bivariate  exponential  (ACBVE)  distri- 
bution used  as  a prior  in  the  exponential  failure  case  may  also  be  used 
here. 

Contours  of  some  posterior  distributions  for  (C,P)  , together 
with  their  modes,  are  shown  in  Figure  3.3(i)  - (vi) . The  contours  la- 
beled A,  B,  C are  the  approximate  50,  75,  and  95%  H.P.D.  regions,  respec- 
tively, obtained  as  described  in  Section  2.1  of  this  chapter.  The  poste- 
rior contours  shown  in  Figure  3.3(i)  - ( i i ) were  obtained  using  a uniform 
prior.  The  posterior  contours  shown  in  Figure  3.3(iii)  - (vi)  were  ob- 
tained using  ACBVE  prior  distributions  with  parameters  given  in  Table  3.4. 
For  each  case  shown  in  Figure  3.3(i)  - (vi)  the  values  a = 0.031  , 
y = 11.375  , A = 0.002  , A2  = 0.25  , M = 47  , and  N = 49  were 

used.  The  two  values  of  b used  to  obtain  the  posterior  contours  shown 
in  Figure  3 . 3 ( i)  - (vi)  are  b = 1.0  , which  is  equivalent  to  assuming 
that  the  failure  times  follow  an  exponential  distribution,  and  b = 0.808 
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TABLE  3.3 

OBSERVED  TIMES  TO  FAILURE  OF  AN  INSULATING  FLUID  AT 
DIFFERENT  VOLTAGE  STRESS  LEVELS  (FROM  NELSON  [15]) 


0.27  0.19  0.35  0.09 

0. AO  0.78  0.59  0.39 

0.69  0.96  0.96  0.A7 

0.79  1.31  0.99  0.73 

2.75  2.78  1.69  0.74 

3.91  3.16  1.97  1.13 

9.88  4.15  2.07  1.40 
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215.10 

32.52 

33.91 

36.71 
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.000  ii)  b = 1.000 


Figure  3.3(i)  - (vi) — Contours  of  the  posterior  distributions 
of  (C,P)  obtained  using  different  priors  with  the  same  value  of  b : Nelson  l 15]  data 
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is  the  value  of  the  least  squares  estimate  obtained  by  Nelson  [15]. 

The  shape  and  location  of  the  contours  of  the  posterior  distributions 
shown  in  Figure  3.3  suggest,  for  the  cases  considered,  that  the  choice 
of  b has  a greater  effect  on  the  posterior  than  the  choice  of 
priors. 

The  Bayesian  framework  has  the  advantage  of  allowing  for  the  ad- 
justment of  the  posterior  in  terms  of  the  shape  parameter  b . That  is, 
by  conditioning  p(C,P/f),b)  on  different  values  for  b , we  may  obtain 
an  indication  of  how  sensitive  inferences  concerning  (C,P)  are  to 
changes  in  b . Three  contours  of  posterior  distributions  of  (C,P)  , 
together  with  their  modes,  for  different  choices  of  b , are  shown  in 
Figure  3.4.  In  each  case,  uniform  priors  were  used  with  a = 0.031  , 
y = 11.  375  , A = 0.002  , A,,  = 0. 25  , M = 47  , and  N = 49  . The 

choice  of  b = 1.0  and  b = 0.808  was  described  above,  b = 0.7709 
and  b = 0.6962  were  obtained  using  the  maximum  likelihood  estimation 
procedures  of  Singpurwalla  and  Al-Khayyal  [21],  b = 0.925  and  0.610 
are  roughly  the  95%  confidence  limits  obtained  by  Nelson  [12],  and 
b = 0. 5 was  arbitrarily  chosen. 

The  contours  shown  in  Figure  3.4  give  an  indication  of  how  in- 
ferences about  (C, P)  are  affected  by  the  choice  of  b . As  b 
changes  within  a small  range,  the  shape  of  the  posterior  contours 
change  significantly.  An  assumption  of  exponentiality  (b  = 1.0)  pro- 
duces a rather  different  result  than  is  obtained  using  the  least  squares 
estimate  (b  = 0.808)  or  the  maximum  likelihood  estimates  (b  = 0.7709 
or  0.6962).  In  particular,  the  dispersion  of  the  joint  posterior  dis- 
tribution of  (C,P)  increases  markedly  as  b decreases  from  1.0  to  0.5. 

Changes  in  the  joint  posterior  mode  for  the  different  values  of 
b are  shown  in  Table  3.5.  With  the  exception  of  b * 0.5  , different 

A 

values  for  C account  for  changes  in  the  mode  of  the  joint  posterior. 

The  effect  of  changes  in  b on  the  distributions  of  C and  P 
may  be  examined  by  obtaining  the  marginal  posterior  distributions  con- 
ditioned on  b : 


r-'ivt 


p(c  /n,b)  = 

m ~ 

and 

p(p  /n,b)  = 

n - 

Plots  of  marginal  distributions  of  C and  P are  shown  in  Figures 
3.5  and  3.6.  These  were  obtained  from  the  joint  posterior  distribu- 
tions with  the.  contours  shown  in  Figure  3.4.  The  marginal  distribu- 
tions also  show  a marked  difference  for  the  different  values  of  b , 
although  the  effect  on  C is  more  pronounced.  The  marginal  distri- 
bution of  C changes  both  dispersion  and  location  as  b decreases 
from  1.0  to  0.5.  The  dispersion  of  the  marginal  distribution  of  P 
changes  significantly  as  b decreases  from  1.0  to  0.5,  but  its  loca- 
tion remains  essentially  the  same. 

We  again  emphasize  in  this  case  that  the  results  obtained  here  ap- 
ply only  for  the  particular  set  of  data  analyzed.  These  results  indicate 
that  inferences  regarding  (C,P)  are  sensitive  (i.e.,  non-robust)  to  the 
value  of  the  Weibull  shape  parameter  b . This  may  be  explained  by  the 
fact  that  the  estimators  have  gamma  distributions  only  when  the  cor- 

rect value  of  b is  known.  The  likelihood  used  in  this  case  is  ob- 
tained by  taking  the  joint  distribution  of  (o.  ,fj  , . . . ,f|  ) conditional 

1 ^ K 

on  a value  for  b . If  the  value  selected  for  b is  not  close  to  the 
correct  value,  we  would  reasonably  expect  inferences  based  on  the  gamma 
assumption  to  be  affected.  By  conditioning  the  posterior  distributions 
on  different  values  of  b , we  obtain  an  indication  of  the  extent  of 
these  effects. 


- 44  - 
N 

) p(C  ,P  /n,b)  , m=0,l,...,M 

i=0  m n ' 

M 

£ p(c  ,P  /h,h)  , n=0, 1 , . . . ,N  . 
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3.2  Prediction  at  Use  Stress 


Obtaining  a distribution  of  predicted  values  of  0 in  the 

Weibull  case  follows  the  same  lines  as  the  exponential  case  discussed 

in  Section  2.2.  That  is,  we  assume  that  a predicted  value  for  9 is 

u 

f 

obtained  as 
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1-1 


0 


u(m,n) 


C (V  /V) 
tn  ii 


where  V is  use  condition  stress  and  V is  given  by  Iquation  0.2.1). 
u 

The  probability  associated  with  0 . . is  then 

r 7 u(m,n) 


p(C  ,P  /f|,b)  . 
tn  n ~ 


The  process  of  obtaining  a distribution  of  predicted  values  and  approxi- 
mate 95%  probability  limits  for  0^  based  on  pCC^.P^/h.b)  is  the  same 

as  the  exponential  case  described  in  Section  2.2. 

Some  distributions  of  predicted  values  for  0 at  V = 20  are 

u u 

shown  in  Figure  3.7.  These  distributions  were  obtained  from  the  pos- 
terior distributions  with  the  contours  shown  in  Figure  3.4.  In  Table 
3.5  are  shown  the  predicted  values  0^  obtained  using  the  modal  values 

(C,P)  and  approximate  95%  probability  limits  for  6^  obtained  from 

the  posterior  distributions  shown  in  Figure  3.4.  The  dispersions  of 
these  prediction  distributions,  as  indicated  by  the  probability  limits, 
are  quite  sensitive  to  the  choice  of  b . However,  the  locations  and 
shapes  of  the  prediction  distributions  are  somewhat  similar. 
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CHAPTER  IV 

SlfMMARY  AND  CONCI.USIONS 


In  this  dissertation  we  have  considered  some  problems  of  statis- 
tical inference  in  accelerated  life  testing.  Bayesian  and  non-Bayesian 
approaches  to  these  problems  were  presented.  In  this  chapter  the  re- 
sults are  summarized  and  conclusions  drawn  concerning  their 
applicability . 

In  the  first  part  of  this  dissertation  we  obtained  orthogonal 
least  squares  estimators  for  the  parameters  of  three  models  used  to 
re-parameter ize  the  scale  parameter  of  an  exponential  distribution:  the 

power  rule  model,  the  Arrhenius  model,  and  the  Eyring  model.  We  estab- 
lished that  these  estimators  have  some  convenient  and,  from  the  classi- 
cal or  non-Bayesian  viewpoint,  desirable  properties. 

(i)  They  are  best  linear  unbiased  and  orthogonal  for  any 
sample  size,  i.e.,  any  number  of  stress  levels,  greater  than  the  number 
of  parameters  being  estimated. 

(ii)  The  exact  variances  and  higher  moments  of  the  estimators 
may  be  readily  obtained  for  any  test  design. 

(iii)  The  least  squares  estimators  are  asymptotically  normally 
distributed.  Further,  the  estimators  may  be  considered  approximately 
normally  distributed  for  tests  involving  small  numbers  of  stress  levels 
provided  a moderate  number  of  failures  is  observed  at  each  stress. 

(iv)  The  least  squares  approach  avoids  the  solution  of  a num- 
ber of  non-linear  equations  which  is  required  by  maximum  likelihood 
estimation.  The  number  of  non-linear  equations  to  be  solved  is  equal 

to  the  number  of  parameters  being  estimated  and  this  could  be  cumbersome. 


49  - 


i 


,y 


A1  though  the  least  squares  estimators  are  not  asymptotically  efficient, 
their  variances  are  exact  anti  only  slightly  greater  than  the  Cramer- 
Rao  lower  bound  for  any  sample  size.  Also,  the  least  squares  estima- 
tors have  zero  covariance  for  any  sample  size  while  the  covariance  of 
the  maximum  likelihood  estimates  is  zero  asymptotically. 

(v)  Best  linear  unbiased  estimators  of  linear  functions  of 
the  unknown  parameters,  such  as  the  log  of  the  exponential  scale  param- 
eter at  use  or  normal  stress,  are  obtained  by  substituting  the  least 
squares  estimators  in  the  function. 

In  the  second  part  of  this  dissertation  we  considered  a Bayesian 
approach  for  the  analysis  of  the  power  rule  model.  Two  cases  were 
studied:  the  power  rule  used  to  re-parameterize  the  scale  parameter  of 

an  exponential  distribution  and  the  scale  parameter  of  a Weibull  distri- 
bution. In  the  Bayesian  analysis  we  obtained  joint  posterior  distribu- 
tions for  the  two  parameters  of  the  power  rule.  The  contours  of  the 
posterior  distributions  shown  in  this  dissertation  indicate  regions  of 
the  parameter  space  which  have  high  posterior  probability  of  containing 
the  true  parameter  values.  The  posterior  distributions  were  obtained 
using  a numerical  procedure  which  involved  a discretization  of  the  prior 
distribution.  In  both  the  exponential  and  Weibull  cases,  posterior  dis- 
tributions were  obtained  using  uniform  and  absolutely  continuous  bivari- 
ate exponential  (ACBVE)  prior  distributions.  Location  parameters  were 
introduced  in  the  ACBVE  distribution  in  order  to  have  more  flexibility 
in  specifying  the  location  of  the  prior  distribution  with  regard  to  the 
region  of  the  parameter  space  for  which  the  likelihood  is  appreciably 
non-zero. 

For  the  exponential  case,  data  from  simulated  accelerated  life 
tests  were  analyzed  with  a uniform  prior  distribution  and  with  ACBVE 
prior  distributions  with  different  parameter  values.  The  contours  of 
the  posterior  distributions  of  the  power  rule  parameters  were  quite  sim- 
ilar due  to  the  dominance  of  the  realized  likelihood  in  a relatively 
small  region  of  the  parameter  space  and  the  lack  of  variation  of  the 
ACBVE  priors  in  this  region.  For  the  Weibull  case,  data  from  acceler- 
ated life  tests  of  an  insulating  fluid  under  voltage  stresses  were 
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analyzed.  Posterior  distributions  for  the  power  rule  parameters  were 
obtained  conditional  on  different  values  for  the  Weibull  shape  param- 
eter using  uniform  and  different  ACBVK  prior  distributions.  The 
Bayesian  approach  has  the  advantage  of  clearly  showing  the  sensitivity 
of  inference  procedures  to  assumptions  about  the  shape  parameter  and 
prior  knowledge.  Different  values  of  the  shape  parameter,  within  a 
narrow  range,  were  shown  to  have  a marked  effect  on  the  posterior  dis- 
tribution of  the  power  rule  parameters.  That  is,  using  the  same  prior 
distribution  for  the  power  rule  parameters,  slight  changes  in  the  shape 
parameter  were  shown  to  cause  significant  changes  in  the  dispersion  and 
location  of  the  posterior  distribution.  Thus,  for  the  data  analyzed, 
we  conclude  that  inferences  about  the  power  rule  parameters  are  not 
robust  with  respect  to  assumptions  about  the  Weibull  shape  parameter. 

The  joint  posterior  distribution  of  the  power  rule  parameters  may 
be  used  to  obtain  inferences  about  functions  of  the  parameters.  The 
posterior  distributions  obtained  in  tli is  dissertation  were  used  to  con- 
struct distributions  of  predicted  values  of  the  exponential  and  Weibull 
scale  parameters  at  use  condition  or  normal  stress.  In  the  exponential 
case  the  distributions  of  predicted  values  for  the  scale  parameter  were 
quite  similar  since  the  posterior  distributions  of  the  power  rule  param- 
eters were  similar.  In  the  Weibull  case  the  distributions  of  predicted 
values  for  the  scale  parameter  reflected  the  sensitivity  to  changes  in 
the  Weibull  shape  parameter  of  the  distribution  of  the  power  rule 
parameters . 

The  question  naturally  arises  as  to  which  of  the  approaches 
presented  here,  the  classical  least  squares  or  the  Bayesian,  should  be 
used  in  practice.  An  answer  to  this  question  would  require  that  basic 
philosophical  issues  regarding  the  differences  between  Bayesian  and  non- 
Bayesian  methods  be  resolved.  In  the  absence  of  a generally  accepted 
resolution  of  these  issues  the  choice  of  what  methods  to  use  depends  on 
the  preferences  of  the  analyst.  A non-Bayesian  in  search  of  a classical 
solution  to  the  problems  considered  here  might  prefer  the  least  squares 
approach  since  it  results  in  estimators  which  have,  from  the  classical 
viewpoint,  desirable  properties  and  presents  no  computational  problems. 


On  the  other  hand,  a Bayesian  would  reject  the  least  squares 
approach,  or  any  other  inference  procedure  that  does  not  result  from  a 
likelihood  function  and  a set  of  prior  distributions,  as  having  defi- 
ciencies in  general  cases  that  may  be  verified  from  an  objective  view- 
point. Indeed,  general  acceptance  of  this  idea  would  resolve  the  ques- 
tion in  favor  of  the  Bayesian  approach  for  all  inference  problems.  At 
the  least,  the  Bayesian  approach  has  some  advantages  that  should  be 
considered : 

(i)  Most  analysts,  however  neutral  philosophically,  have  some 
prior  knowledge  or  beliefs  about  their  problems.  The  Bayesian  framework 
allows  for  the  quantitative  expression  of  this  prior  knowledge. 

(ii)  The  model  or  likelihood  function  through  which  the  data 
act  to  modify  prior  knowledge  must  be  specified  explicitly.  This  re- 
quires that  careful  consideration  be  given  to  the  process  generating  the 
data  while  also  facilitating  the  exploration  of  the  consequences  of  al- 
ternative assumptions  regarding  the  process. 

(iii)  In  the  Bayesian  framework  inferences  are  expressed  in 
terms  of  posterior  distributions.  This  gives  a probability  interpreta- 
tion to  the  inferences  that  is  not  possible  using  the  classical  confidence 
intervals. 
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AITKNIHX  A 

SHAPE  FACTORS  OF  THE  LEAST  SQUARES  ESTIMATES  OF 
THE  POWER  RULE  AND  ARRHENIUS  MODEL  PARAMETERS 

In  this  appendix,  we  derive  expressions  for  the  skewness  and 
kurtosis  of  the  least  squares  estimates  Bq  and  for  both  the 

power  rule  and  the  Arrhenius  model  parameters.  The  skewness  and  kurto- 
sis, or  shape  factors,  give  an  indication  of  the  extent  to  which  the 
distributions  of  the  estimates  deviate  from  normality.  Also,  since  the 
distributions  of  Bq  and  B^  approach  the  normal  as  k , the  number 

of  stress  levels,  increases,  the  shape  factors  give  an  indication  of 
how  rapid  the  approach  is. 

The  least  squares  estimates  for  the  power  rule  and  the  Arrhenius 
* /\  - , 
models,  B„  and  B,  , are  linear  combinations  of  the  observations  iz.l 
0 1 1 

which  are  k independent  non- identically  distributed  log  gamma  random 

3/2 

variables  with  Var(z.)  = K„  , skewness  y,  = K„/K~  , and  kurtosis 

i 2 1 , z . 3 2 

Y2,z  = VK2  * where 

1 

Ks+1  = ^(S)(r1)  ’ s > 0 

is  the  (s+l)st  cumulant  of  a log  gamma  random  variable  with  shape  param- 

( s ) 

eter  r^  , and  ip  (•)  is  the  sth  derivative  of  the  digamma  function. 

Thus,  using  formulas  given  by  Scheffe  ([18],  page  332),  we  may  compute 
the  skewness  and  kurtosis  of  B^  and  B^  or  any  other  linear  combina- 
tion of  the  {z^}  • 

k 

In  general,  for  a linear  combination  £ = £ c.y.  °f  k inde- 

i=l  1 1 

2 

pendent  random  variables  {y^}  with  variance  , skewness  y^  ^ and 
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kurtosis  y . , the  skewness  and  kurtosis  of  £ are 

Z , 1 


V 3/2 

'l,S  * <1>ai  h,i 


where  (±)  is  the  sign  of  , and 


» r = l V2  i ’ 
i=1  1 


2 2/  22 

where  a.  = c.a.  / J*  c.o.  . The  a.  are  the  proportion  of  the  vari- 

1 1 1 / 1 1 1 rr 

1 k 

ance  of  E,  contributed  by  the  ith  term,  and  thus  £ a.  = 1 , For 

i=l  1 

§0  ’ ai,30  = WiVi=l  W_il  and  f0r  §1  ’ = Wilxi  /ill  Wilxi  ’ 

where  the  w^  are  given  by  Equation  (2.2.4)  and  x.  is  given  by  Equa- 
tion (2.2.7)  for  the  power  rule  and  Equation  (2.3.3)  for  the  Arrhenius 
model . 

We  may  note  that  the  effects  of  observing  a small  number  of  fail- 
ures at  one  or  more  stress  levels  in  an  accelerated  test  may  be  offset 
by  observing  a larger  number  of  failures  at  other  stress  levels.  That 
is,  since  the  shape  factors  of  the  estimators  are  weighted  averages  of 
the  shape  factors  of  the  individual  observations,  the  effect  of  obtain- 
ing a few  failures  at  one  or  more  stress  levels  may  be  balanced  by  ob- 
serving larger  numbers  of  failures  at  other  stress  levels. 

The  following  results  may  be  shown  to  hold  if  the  r^  , the  num- 
ber of  failures  observed  at  the  ith  stress  level,  are  unequal.  We  will 
assume,  however,  that  r^  = r for  all  i since  this  results  in  less 
complicated  expressions. 


When  r = r , i=l,..^,k  , the  variance,  skewness  and  kurtosis  of 

the  z.'s  are  equal  and,  using  the  notation  w.  = w , y.  = y,  and 
l l i,z  1 

> \ 

y»  - y~  , i=l,...,k  , it  follows  that  ot.  g = 1/k  , i=l,...,k  , 

Z,Zi  1 1,fS0 


| r-  -«•» 


• I V . 


1 
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and  y 


-1 


1,B„ 


= k , and  Y2  £ = k y2  • Clearly,  as 


k Is  increased  with  r fixed,  the  shape  factors  for  fig  approach 
the  normal  distribution  value  of  zero.  Using  the  approximations  given 


by  Johnson  and  Kotz  ( [8]  f page  196)  for  the  shape  factors  of  a log 

1 

i | 2 

gamma  random  variable  with  shape  parameter 


r ’ Yi  * - (r  - i) 


and 


Y2  * 2 


(r  - |)  , it  i« 


is  clear  that  if  k is  held  fixed,  increasing 
r also  causes  the  shape  factors  to  approach  zero. 


For  B.  , if  r.  = r , i=l,...,k  , a " 

1 1 1»“1 


2 2 

x . / £ x . , i=l , . . . ,k  , 

Xl  j=l  1 


y / 2 / y 2V/2  . y / 2 / ^ 2V 

= 1 I x ■ / s x • I Y.  and  = £ I x . / £ x.  Y,  • 

1'h  i=i\  x/  j=i  V 1 2’6i  i=i\  x/ j=i  J/  2 


so  that  y 


Again,  using  the  approximations  for  y^  and  y2  , it  is  clear  that 
holding  k fixed,  increasing  r causes  the  shape  factors  for  8^  to 
approach  zero.  Holding  r fixed  as  k -*■  °°  , the  rate  of  approach  to 


zero  of  the  shape  factors  for  8^  will  be  similar  to  that  for  6q  if 


the  a . 


1,6. 


are  approximately  equal  for  all  i . In  many  cases  of  the 


power  rule  and  the  Arrhenius  models,  however,  one  or  more  of  the  x^ 
has  a limiting  value  which  is  disproportionately  larger  than  the  others 


and  thus  k must  be  quite  large  before  the  corresponding  a.  " 

i,6. 


becomes  small,  which  in  turn  allows 


r J'2 

1 a-i  ft 

i=l  1,B1 


and 


E a 


i.Bi 


to 


i=l  “"l 

approach  zero.  Although  in  practice  k is  usually  rather  small,  the 


shape  factors  for  both  Bq  and  B^  can  be  made  reasonably  close  to 

the  normal  distribution  values  for  moderate  size  r , as  shown  in 
Table  A. 1 . 


i 
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TABLE  A. 1 


SHAPE  FACTORS  FOR  THE  ORTHOGONAL  LEAST  SQUARES  ESTIMATORS 
OF  THE  POWER  RULE  ANI)  ARRHENIUS  PARAMETERS 


k 

= 5 

k 

= 10 

r 

b,i0 

^2  S 
1 o 

v 

'-60 

r2,e0 

10 

-.14510 

.04211 

-.10260 

.02105 

15 

-.11744 

.02759 

-.08305 

.01379 

20 

-.10127 

.02051 

-.07161 

.01026 

r 

Y1 

’P1 

Y2 

’ei 

10 

- . 32444S 

. 21053K 

15 

- . 26261S 

. 13793K 

20 

-.22646 S 

. 10256K. 

where 


r 

k 

S 

K 


number  of  failures  observed 
at  each  stress  level 


number  of  stress  levels 
£ 3/2 

Z 01  i ft 

i=l  ,ei 

? 2 . 

i=l  1,B1 
2 < k 2 

x.  I E x , i=l, . . . ,k  , 
! j=i  3 


where  the  are  defined  by  Equation  (2.2.7)  for  the  power  rule 

and  Equation  (2.3.3)  for  the  Arrhenius  model, 
k 

Note:  Since  I a . ' = 1 , S,K  < 1 . 

i=l  1,1 1 
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APPENDIX  B 

SHAPE  FACTORS  OF  THE  LEAST  SQUARES  ESTIMATES 
OF  THE  EYRING  MODEL  PARAMETERS 


In  this  appendix,  we  derive  expressions  for  the  skewness  and  kur- 
tosis  of  the  orthogonal  least  squares  estimates  of  the  Eyring  parameters, 
S'  = (0^ , S>2 » 3^)  given  by  Equation  (2.6.8).  The  skewness  and  kurto- 

sis,  or  shape  factors,  give  an  indication  of  the  extent  to  which  the  dis- 
tributions of  the  estimates  deviate  from  normality.  Also,  since  the  dis- 
tributions of  the  estimates  approach  the  normal  as  k , the  number  of 
stress  levels  increases,  the  shape  factors  give  an  indication  of  how 
rapid  the  approach  is. 


The  orthogonal  least  squares  estimates  of  the  Eyring  model  param- 
eters are  linear  combinations  of  the  observations  (z which  are  k 


independent  non-identically  distributed  log  gamma  random  variables  with 

3/2  2 

Var(z.)  = K.  , skewness  v = K./K„  and  kurtosis  y„  = K, /K„  , 
l 2 1 , z . 3 2 2 , z . 4 2 


where 


K 


s+1 


0 


s t 

is  the  (s+1)  cumulant  of  a log  gamma  random  variable  with  shape  pa- 

rameter r^  and  i]/S^(*)  is  the  s^  derivative  of  the  digamma  func- 
tion. Thus,  using  formulas  given  by  Sheffe  ([18],  p.  332),  we  may  com- 
pute the  skewness  and  kurtosis  of  the  least  squares  estimates  or  any 
other  linear  combination  of  the  {z^}  • 


The  general  formulas  for  the  shape  factors  of  a linear  combination 
of  independent  random  variables  are  given  in  Appendix  A.  For  the  orthog- 
onal least  squares  estimates  of  the  Eyring  model  parameters 
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3/2 

- £ ai,3.  Y1,Z. 


, 3=1, • ••,4 


2 

y n n ~ Olj  n y p „ 

2,3j  i=l  ,3j  2’  i 


, .1=1,2, 3,4 


where 


/ k 

wTL  / E wT1 
1 /j=l  J 


, 1=1, 


,k 


„ „ 2 —1  / 2,  2-1 

x’32  12  " / j=l  j2  1 


,k 


2 -1  it  2-1  . . 

x = x w.  / E x w.  , i=l,. 

i,e3  13  1 / j=1  ,2  j 


.,k 


2 -1  / _ 2-1 

t.  . = x.,w.  / £ x..w.  , i=l,. 

i,64  i / j=i  J 


, k 


w..  = 


, i 1 , . • . , k 


and  x^0  , x^  , x^  for  1=1,..., k are  given  by  Equations  (2.6.9), 
(2.6.10),  and  (2.6.11),  respectively. 


The  following  results,  which  relate  to  the  behavior  of  the 
shape  factors  of  the  estimates  as  the  number  of  stress  levels  and  the 
number  of  failures  observed  at  each  stress  increase,  may  be  shown  to 
hold  if  the  r.  , the  number  of  failures  observed  at  the  i1"*1  stress 

l 

level,  are  unequal.  However,  assuming  r^  = r for  all  i results 
in  less  complicated  expressions  and  causes  no  loss  of  generality. 


When  r^  = r , i=l,...,k  , the  variance,  skewness  and  kurtosis 

of  the  z.'s  are  equal.  Using  the  notation  w.  = w , y,  = Y, 

i l 1 ,z  1 

and  z = expressions  for  the  skewness  and  kurtosis  of  the 

~2  , -1 

estimates  are  obtained  as  y,  ' = k y.  , y»  ' = k y_  and 

mPi  J-  *-  » Pi  ^ 
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'i-s. = Ti  a fall **■) 

‘ Y2  ‘ml  *L) 


for  m = 2,3,4 


Clearly,  as  k is  increased  with  r held  fixed,  the  shape  factors  for 
3^  approach  the  normal  distribution  values  of  zero.  Using  the  approxi- 
mations given  by  Johnson  and  Kotz  ([8],  p.  196)  for  the  shape  factors 
of  a log  gamma  random  variable  with  shape  parameter  r , 

1 - — 1-1 

W -(r  - — ) 2 and  Y2~2(r  - — ) , it  is  clear  that  if  k is  held 

fixed,  increasing  r also  causes  the  shape  factors  of  8^  to  approach 

zero.  For  3 , m = 2,3,4  , again  using  the  approximations  for  y,  and 

m 1 

>2  , it  is  clear  that  increasing  r while  holding  k fixed  causes 

y ~ and  y ~ to  approach  zero.  Holding  r fixed  as  k -*•<*>  , 

1,8  Mi,, 

m m 

the  rate  of  approach  to  zero  of  the  shape  factors  for  8 , m = 2,3,4 

m 

will  be  similar  to  that  for  the  shape  factors  of  8^  if  the 

2 / k 2 

Y . 2 ~ x.  I T.  x.  , j=l,...,k  , are  approximately  equal  for  all  j . 

J’6m  Jm/i=l  lm 

2 

In  many  cases  of  the  Eyring  model,  however,  one  or  more  of  the  x, 

has  a limiting  value  which  is  disproportionately  larger  than  the  others 
and  thus  k must  be  quite  large  before  the  corresponding  a.  ' 

3 * 

k k m 

3/2  2 

becomes  small,  which  in  turn  allows  E a . ~ and  £ a.  s to 

j=l  J’Pm  J-l  J’6m 

approach  zero.  Although  in  practice  k is  usually  rather  small,  the 

shape  factors  for  8 , m = 1,2, 3, 4 , can  be  made  reasonably  close  to 

m 

the  normal  distribution  values  for  moderate  size  r , as  shown  in 
Table  B.l. 


» ► . 


r = number  of  failures  observed  at 
each  stress  level 

k = number  of  stress  levels 


s 

m 

k 

= Z 

j=l 

a3/? 

J’em 

, m = 2,3,4 

U 

m 

k 

= Z 

3=1 

2 A 

a . ^ 

m 

, m = 2,3,4 

2 

= x . 
3, 

/ f 2 

/ E x . 
m / . , im 

/ 1=1 

, 3 = l,...,k 

m = 2,3,4 

With  X.  ry  , 

J ’ ^ 

X. 

3.3 

• X3,4  ’ j 

= 1, . . . ,k  given  by 

Equations  (2.6.9),  (2.6.10)  and  (2.6.11),  respectively. 


k 

Note:  Z a'  =1,  S,U  <1  for  m=  2,3,4 

1 , p mm 
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